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Abstract 

A review of electronic dynamics of single-impurity and many-impurity Anderson models is 
contained in this report. Those models are used widely for many of the applications in diverse 
fields of interest, such as surface physics, theory of chemisorption and adsorbate reactions 
on metal surfaces, physics of intermediate valence systems, theory of heavy fermions, physics 
of quantum dots and other nanostructures. While standard treatments are generally based 
on perturbation methods, our approach is based on the non-perturbative technique for the 
thermodynamic Green functions. The method of the irreducible Green functions is used as the 
basic tool. This irreducible Green functions method allows one to describe the quasiparticle 
spectra with damping of the strongly correlated electron systems in a very general and natural 
way and to construct the relevant dynamical solution in a self-consistent way on the level 
of Dyson equation without decoupling the chain of the equations of motion for the Green 
functions. The subject matter includes the improved interpolating solution of the Anderson 
model. It was shown that an interpolating approximation, which simultaneously reproduces the 
weak-coupling limit up to second order in the interaction strength U and the strong-coupling 
limit up to second order in the hybridization V (and thus also fulfils the atomic limit) can be 
formulated self-consistently. This approach offers a new way for the systematic construction 
of approximate interpolation dynamical solutions of strongly correlated electron systems. 


Keywords: Many-particle interacting systems; statistical physics; physical chemistry; surface 
physics; theory of chemisorption; fundamental aspects of catalysis; hybridizing localized and 
itinerant electrons; the single-impurity, two-impurity and periodic Anderson model; the elec¬ 
tronic quasiparticle dynamics; the Green functions approach; Dyson equation; quasiparticle 
damping. 
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1 Introduction 


In this review we discuss the many-body quasiparticle dynamics of the Anderson impurity model 1 2 
and its generalization^® in the framework of the equation-of-motion method®EM 7][8l[9][To][TT][T2] 
hnite temperatures. 

The studies of strongly correlated electrons in solids and their quasiparticle dynamics are inten¬ 
sively explored subjects in solid state physicsEE3E II Electronic dynamics in the bulk and at 
the surface of solid materials are well known to play a key role in a variety of physical and chem¬ 
ical phenomena!®®^ 617 18 19 20 One of the main aspects of such studies is the interaction of 
low-energy electrons with solids, where the calculations of inelastic lifetimes of both low-energy 
electrons in bulk materials and image-potential states at metal surfaces are highly actual problems. 
The calculations of inelastic lifetimes was made as a rule in a model of the homogeneous electron 
gasf 1 ' 5 by using various approximate representations of the electronic response of the medium. 
Band-structure calculations, which have been carried out in literature may give a partial informa¬ 
tion only. 

The band-structure approacrf®®®® suffers from well known limitations. It cannot be validated 
in full measure in the case of very narrow bands and strongly correlated localized electrons. 5 24 25 
An alternative approach is connected with using correlated fermion lattice models, like Ander- 

The principal importance of this approach is related with the dual character of electrons in dilute 
magnetic alloys EIMM33I3I11M11401 j n transition metal oxides EM3HII21H31H31H51ISll3ZlL48j|49j[So] intermediate- 
valence 


heavy fermionsP l &1 || 54 || 55 |l &6 1 high-T c superconductors, 9 57 - etc. In these ma¬ 
terials electrons exhibit both localized and delocalized features. 5 41 For example in paper 5 ® the 
electronic structure in solid phases of plutonium was discussed. The electrons in the outermost 
orbitals of plutonium show qualities of both atomic and metallic electrons. The metallic aspects 
of electrons and the electron duality that effect the electronic, magnetic and other properties of 
elements were manifested clearly. 

The basic models to describe correlated electron systems are the single-impurity Anderson model 
(SIAM),®®®®® 8 59 p er iodic Anderson model (PAM) 60 and the Hubbard model which exhibit 
the key physical feature, i.e., the competition between kinetic energy (itinerant) and potential 
energy (localized) effects!®®®] 

Indeed, the standard approach which is valid mainly for the simple and noble metals is provided 
by the band theory formalism for the calculation of the electronic structure of solids. For a better 
understanding of how structure and properties of solids may be related the chemically insightful 
concept of orbital interaction and the essential machinery of band theory should be taken into 
account 2 ® to reveal links between the crystal and electronic structure of periodic systems. In such 
a way, it was possible shown 25 how important tools for understanding properties of solids like the 
density of states, the Fermi surface etc., can be qualitatively formulated and used to rationalize 
experimental observations. It was shown that extensive use of the orbital interaction approach 
appears to be a very efficient way of building bridges between physically and chemically based 
notions to understand the structure and properties of solids. 

The Anderson and Hubbard models found a lot of applications in studies of surface physics,®® 
theory of chemisorption and adsorption!®®®®® and various aspects of physics of quantum 
dots.® 6 ' 67 68 69 70 71 < 2 , 73,74 75 . 76,77 However in spite of many theoretical efforts a fully satisfac¬ 
tory solution of the dynamical problem is still missing. The Bethe-ansatz solution of the SIAM 
allows for the determination of the ground state and thermodynamic static properties, but it does 
not allow for a determination of the dynamical properties. For their understanding the develop¬ 
ment of improved and reliable approximations is still justified and desirable. In this context it 
is of interest to consider an interpolating and improved interpolating approximations which were 
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proposed in the papers!* 3,7,8 9,10,11 12 We will show that a self-consistent approximation for the 
SIAM can be formulated which reproduces all relevant exactly solvable limits and interpolates 
between the strong- and weak-coupling limit. 

In connection with the dynamical properties the one-particle Green function is the basic quan¬ 
tity to be calculated. Subject of this survey is primarily devoted to the analysis of the relevant 
many-body dynamic solution of the single impurity Anderson model and its correct functional 
structure. We wish to emphasize that the correct functional structure actually arises both from 
the self-consistent many-body approach and intrinsic nature of the model itself. The important 
representative quantity is the spectral intensity of the Green function at low energy and low tem¬ 
perature. Hence, it is desirable to have a consistent and closed analytic representation for the 
one-particle Green function of SIAM. The paperJ 6,7 ®®!^^] dearly show the importance of 
the calculation of the Green function and spectral densities for SIAM and the many-impurity 
Anderson model in a self-consistent way. 

In this terse overview the problem of consistent analytic description of the many-body dynam¬ 
ics of SIAM is analyzed in the framework of the equation-of-motion approach for double-time 
thermodynamic Green functions. 4,5 In addition to the irreducible Green functions (IGF) ap¬ 
proach, 32,78,79 80,81,82 83,84,85 we use a new exact identity EEUE] relating the one-particle and 
many-particle Green functions. Using this identity, it was possible to formulate a consistent and 
general scheme for construction of generalized solutions of the Anderson model. A new approach 
for the complex expansion for the single-particle propagator in terms of Coulomb repulsion U 
and hybridization V is discussed as well. Using the exact identity, an essentially new many-body 
dynamic solution of SIAM was derived. 

2 Hamiltonian of the Models 

2.1 Single-impurity Anderson model (SIAM) 

The Hamiltonian of the SIAM can be written in the form 

H = Y e k C \<j C kc + Y E Oafo a fo<T + (!) 

k<j a 

\ Y fLfoafo-afo-cr + Y V k^lj 0cT + 

a kcr 


where ct and are the creation operators for conduction and localized electrons; er is the 
conduction electron dispersion, Eo a is the localized (f-) electron energy level and U is the intra- 
atomic Coulomb interaction at the impurity site. represents the s — / hybridization. In the 
following consideration we will omit the vector notation for the sake of brevity. 

2.2 Periodic Anderson Model (PAM) 

Let us now consider a lattice generalization of SIAM, the so-called periodic Anderson model 
(PAM). The basic assumption of the periodic impurity Anderson model is the presence of two well- 
defined subsystems, i.e. the Fermi sea of nearly free conduction electrons and the localized impurity 
orbitals embedded into the continuum of conduction electron states (in rare-earth compounds, for 
instance, the continuum is actually a mixture of s, p, and d states, and the localized orbitals are 
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( 2 ) 


/ states). The simplest form of PAM 

* = E tkCfca-Cka ^ ^ -^0 ffia U/% ^ ^ Tl^cr^'z—cr H - 

ka ia ia 

- 7 = y2(exp(ikRi)cl a f i(T + exp(-ikRi)fJ a c ka ) 

ika 

assumes a one-electron energy level Eq, hybridization interaction V, and the Coulomb interaction 
U at each lattice site. Using the transformation 

C L = ~7^ J2 ex PHkRj)c] a ; c ka = -j= exp{ikRj)c ja (3) 

V j V 3 

the Hamiltonian (J3J) can be rewritten in the Wannier representation: 

H = ^2 tij c ia C j° + y! EofJ a fi a + U /2 riiaTli-a + (4) 

ij a ia icy 

4^ + f iijCicr ) ■ 

i<j 

If one retains the fc-dependence of the hybridization matrix element V k in Q, the last term in the 
PAM Hamiltonian describing the hybridization interaction between the localized impurity states 
and extended conduction states and containing the essence of a specificity of the Anderson model, 
is as follows 

*v(4t ficr + ./i c ^); V H = v kexp[ik{Rj - R { )]. (5) 

ija k 

The on-site hybridization Vu is equal to zero for symmetry reasons. Hence the Hamiltonian of 
PAM in the Bloch representation takes the form 

H = ^2 e k c k(T c ka + ^2 E kf ka .fka + U/2 ^ n ia ni- a + (6) 

key icy icy 

y! Vk(c\ a fka + /fco-Cfco-)- 

key 

Note that as compared to the SIAM, the PAM has its own specific features. This can lead to 
peculiar magnetic properties for concentrated rare-earth systems where the criterion for magnetic 
ordering depends on the competition between indirect RKKY-type interaction ! 13 (not included 
into SIAM) and the Kondo-type singlet-site screening (contained in SIAM). 

The inclusion of inter-impurity correlations makes the problem even more difficult. Since these 
inter-impurity effects play an essential role in physical behaviour of real systems , 8 81 it is instructive 
to consider the two-impurity Anderson model (TIAM) too. 

2.3 Two-Impurity Anderson Model (TIAM) 

The two-impurity Anderson model was considered by Alexander and Anderson . 3 They put forward 
a theory which introduces the impurity-impurity interaction within a game of parameters. 

The Hamiltonian of TIAM reads 

H = 'Y] tjjC^Cjg + ^2 E 0ifiafi<r + U / 2 ^2 n io n i-<r + ( 7 ) 

ija 2 = 1 , 2a i=l,2a 

J2( V ^ C lfi* + Vikflcia) + J2^fU^ + V ^flha) 

ia a 
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where E o* are the position energies of localized states (for simplicity, we consider identical impuri¬ 
ties and s-type i.e. non-degenerate) orbitals: Eq\ = £02 = £o- The hybridization matrix element 
Vik was discussed in in detail in Ref . 8 As for the TIAM, the situation with the right definition of 
the parameters V \2 and Vik is not very clear. The definition of V \2 in 3 is the following: 

V 12 = 4 = J HfM^dr . (8) 

Note that Hf is without ”H-F” (Hartree-Fock) mark. The essentially local character of the 
Hamiltonian Hf clearly shows that V \2 describes the direct coupling between nearest neighboring 
sites (for a detailed discussion see Ref . 8 where the hierarchy of the Anderson models was discussed 
too). 

3 The Method of Irreducible Green Functions 

When working with infinite hierarchies of equations for Green functions 4 3 the main problem is 
finding the methods for their efficient decoupling, with the aim of obtaining a closed system of 
equations, which determine the Green functions. In the papers 32 78 79 80 81 devoted to studies 
of lattice systems of interacting fermions it was shown that for a wide range of problems in 
statistical mechanics 82 83 84 85 and theory of condensed matter one can outline a fairly systematic 
recipe for constructing approximate solutions in the framework of irreducible Green functions 
method. Within this approach one can look from a unified point of view at the main problems of 
fundamental characters arising in the method of two-time temperature Green functions. 

The method of irreducible Green functions is a useful reformulation of the ordinary Bogoliubov- 
Tyablikov method of equations of motion . 4 5 The constructive idea can be summarized as follows. 
During calculations of single-particle characteristics of the system (the spectrum of quasiparticle 
excitations, the density of states, and others) it is convenient to begin from writing down Green 
function as a formal solution of the Dyson equation. This will allow one to perform the necessary 
decoupling of many-particle correlation functions in the mass operator. This way one can to 
control the decoupling procedure conditionally, by analogy with the diagrammatic approach. In 
this approach the infinite hierarchy of coupled equations for correlation functions is reduced to a 
few relatively simple equations that effectively take into account the essential information on the 
system under consideration, which determine the special features of this concrete problem. 

It is necessary to stress that the structure of solutions obtained in the framework of irreducible 
Green functions method is very sensitive to the order of equations for Green function ^ 4 in which 
irreducible parts are separated. This in turn determines the character of the approximate solutions 


constructed on the basis of the exact representation. 

Now we introduce the retarded, advanced, and causal Green function: 

G r (A, B,t — t') = (( A(t),B(t'))) r = -ie(t-t')([A(t),B(t%), 77 = ±, (9) 

G a (A, B- t-t') = (( A(t ), B(t'))) a = W(t' - t)([A(t),B(t%), 77 = ±, (10) 

G C (A, B; t-t') = (( A(t),B(t'))) c = iT{A{t)B{t')) = (11) 

i6(t — t'){A{t)B{t')) + 7 ji0(t' — t)(B(t')A(t)), 77 = ±. 

Here (...) is the average over the grand canonical ensemble, 8(t) is the Heaviside step function; 
the square brackets denote either commutator or anticommutator (77 = ±): 

[A, B]_ v = AB — rjBA. (12) 


An important ingredient for Green function application is their temporal evolution. In order to 
derive the corresponding evolution’s equation, one has to differentiate Green function over one of 
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its arguments. In order to clarify the above general description, let us consider the equations of 
motion for the retarded Green function of the form 


u>G{ U ) = ([A,Al] v ) + (([A,HUA')) u . 

The irreducible (ir) Green function is defined by 

The unknown constant z is found from the condition 

{[W[A,H]-,A%) = 0. 

It is worth noting that instead of finding the irreducible part of Green function 


(13) 


(14) 


(15) 


(16) 


one can absolutely equivalently consider the irreducible operators 

= (17) 

Therefore, we will use both the notation (W((^4|S))) and (((v4)h r )|5))), whichever is more con¬ 
venient and compact. Equation (1151) implies 


z = 


([[A,H]-,A^]j 


([A,Ai]r 


Mi 

Mo' 


(18) 


Here, Mo and Mi are the zero and first moments of the spectral density.®® Green function is 
called irreducible (i.e. impossible to reduce to a desired, simpler, or smaller form or amount) if 
it cannot be turned into a lower order Green function via decoupling. The well-known objects in 
statistical physics are irreducible correlation functions. In the framework of the diagram technique 
the irreducible vertices are a set of graphs, which cannot be cut along a single line. The definition 
m translates these notions to the language of retarded and advanced Green functions. We 
attribute all the mean-field renormalizations that are separated by Eq. da to Green function 
within a generalized mean field approximation 


GV) = 


([A,A% 

(uj-z) 


(19) 


For calculating Green function (ITU) . ( ir '>(([A, H]_(t), A^(t'))), we make use of differentiation over 
the second time t'. Analogously to Eq. (|14l) we separate the irreducible part from the obtained 
equation and find 

G( w) = G°(w) + G°(u)P(lj)G°(u). (20) 


Here, we introduced the scattering operator 




\-l 


( 21 ) 


In complete analogy with the diagram technique one can use the structure of Eq. (1201) to define 
the mass operator M: 

P = M + MG°P. (22) 
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As a result we obtain the exact Dyson equation (we did not perform any decoupling yet) for 
two-time temperature Green functions: 

G = G° + G°MG. (23) 

According to Eq. (122p . the mass operator M (also known as the self-energy operator) can be 
expressed in terms of the proper (called connected within the diagram technique) part of the 
many-particle irreducible Green function. This operator describes inelastic scattering processes, 
which lead to damping and to additional renormalization of the frequency of self-consistent quasi¬ 
particle excitations. One has to note that there is quite a subtle distinction between the operators 
P and M . Both operators are solutions of two different integral equations given by Eqs. (1221) and 
(1231) . respectively. However, only the Dyson equation ()23[) allows one to write down the following 
formal solution for the Green function: 

G = [{G 0 )- 1 - M]- 1 . (24) 

This fundamental relationship can be considered as an alternative form of the Dyson equation, 
and as the definition of the mass operator under the condition that the Green function within 
the generalized mean-field approximation, G°, was appropriately defined using the equation 

G°G~ 1 + G°M = 1. (25) 

In contrast, the operator P does not satisfy Eq. (1251) . Instead we have 

(G 0 )' 1 - G" 1 = PG°G~ 1 . (26) 

Thus, it is the functional structure of Eq. (124[) that determines the essential differences between 
the operators P and M. To be absolutely precise, the definition (|22l) has a symbolic character. 
It is assumed there that due to the similar structure of equations ([9]) - (fl2l) defining all three 
types of Green functions, one can use the causal Green functions at all stages of calculation, thus 
confirming the sensibility of the definition (1221) . Therefore, one should rather use the phrase ”an 
analogue of the Dyson equation”. Below we will omit this stipulation, because it will not lead to 
misunderstandings. One has to stress that the above definition of irreducible parts of the Green 
function (irreducible operators) is nothing but a general scheme. The specific way of introducing 
the irreducible parts of the Green function depends on the concrete form of the operator A on the 
type of the Hamiltonian, and on the problem under investigation. 

Thus, we managed to reduce the derivation of the complete Green function to calculation of the 
Green function in the generalized mean-field approximation and with the generalized mass oper¬ 
ator. The essential part of the above approach is that the approximate solutions are constructed 
not via decoupling of the equation-of-motion hierarchy, but via choosing the functional form of the 
mass operator in an appropriate self-consistent form. That is, by looking for approximations of 
the form M ~ P[G}. Note that the exact functional structure of the one-particle Green function 
(|24D is preserved in this approach, which is quite an essential advantage in comparison to the 
standard decoupling schemes. 

4 The Irreducible Green Functions Method and SIAM 

After discussing some of the basic facts about the correct functional structure of the relevant 
dynamic solution of correlated electron models we are looking for, described in previous Chapter, 
we give a similar consideration for SIAM. It was shown in Refs. EKEHEI using the minimal algebra 
of relevant operators, that the construction of the generalized mean fields for SIAM is quite 
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nontrivial for the strongly correlated case, and it is rather difficult to get it from an intuitive 
physical point of view. 

To proceed let us consider first the following matrix Green function 


Q( u \ = ( «Cfer|cL» ((cka\foa))\ 
\((/o<rl c L» {(foa\fo„))) ' 


(27) 


Performing the first-time differentiation and defining the irreducible Green function 

~ (^0-ct)((/oct|/o,t))u;) 


(28) 


we obtain the following equation of motion in the matrix form 

J2 Fp(u)G P (u) = i + UD« r \u), (29) 

p 

where all definitions are rather evident. Proceeding further with the IGF technique, the equation 
of motion (129(1 may be rewritten exactly in the form of the Dyson equation 

G[u) = G°{u) + G°{u)M(u)G(uj). (30) 

The generalized mean field Green function G° is defined by 


Y J F p {u)G°p(u) = i. 


(31) 


The explicit solutions for diagonal elements of G° are 

((fo*\fi))Z = (w - E 0a - Un. a - S(u)j) 

(M c L»° = ( w -- u _J 0tT [ Un _ 

where 


-l 


-i 


S(w) = 


u -e k ' 


(32) 

(33) 

(34) 


The mass or self-energy operator, which describes inelastic scattering processes, has the following 
matrix form 

"0 0 


Af( u) = 


0 Moo-y 


where 


ip) 


M 0a = U 2 (^ ir \(f 0a n 0 -Mlrio-*))% r) 

From the formal solution of the Dyson equation (13011 one obtains 

((foa\fl))u, = (w - E 0a - Un-a - M 0a - S( w)) 

„ | t u _ ( _ l^fcl 2 

(( c ka l c far»w u _ e Q(t _ Un-a - Mov 


-1 


-1 


(35) 

(36) 

(37) 

(38) 
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To calculate the self-energy in a self-consistent way, we have to approximate it by lower-order Green 
functions. Let us start by analogy with the Hubbard model with a pair-type approximatioii 32 7981 


[n(w 2 )n(w 3 ) + n(u i)(l 


M 0ff M = U 2 


duj±duj2du}3 

-x 

OJ + 0J\ — 0J 2 — OJ 3 


n(u 2 ) ~ n(uJ3))}g 0 -a(^l)g0a(^2)g0-a(^3), 


(39) 


where we used the notation 

go*(u) = — ^Im ((/Ocrl/ocr))^- (40) 

The equations (1301) and (f39l) constitute a closed self-consistent system of equations for the single¬ 
electron Green function for SIAM model, but only for weakly correlated case. In principle, we can 
use, on the r.h.s. of Eg. (1391) . any workable first iteration-step form of the Green function and find 
a solution by repeated iteration. If we take for the first iteration step the expression 


goa(u) « 6(u - E 0a - Un - a ), 


(41) 


we get, for the self-energy, the explicit expression 

2 n(E 0a + Un-a)( 1 - n(E 0a + Un - a )) 


M 0a (u) = U 


^ E()cr U Tl—fj 

= U 2 Q. a (l - Q_ a )G°(u), 


(42) 


where 

Q-a = n(E 0a + Un- a ), n(E) = {exp[(£ - g)/k B T] + l} -1 . (43) 

This is the well-known atomic limit of the self-energy. 7 10 11 8 

Let us try again another type of the approximation for M. The approximation which we will use 
reflects the interference between the one-particle branch and the collective one 

(fo*(t)f ] 0 _ a (t)fo-a(t)fo-a /<D M ~ 

{foA^fo^ino-ai^no-a} + 

(/o t _ ff (t)/o- ( 7)(/o t ^ (T (^)/o ( 7(t)/o t CT /o- CT ) + 

(fl^t)fo-a)(fo-At)f 0a (t)flfl a ). (44) 


If we retain only the first term in (1441) and make use of the same iteration as in (1411) . we obtain 

(45) 


n/r , X _ TT 2 (! - n ( E 0 c + Un- a )) 1 

Af()<r(Uq ^ U \flQ—(jTlQ— 0 


oj — Eotj — Un- 
If we retain the second term in (|44j) . we obtain 

r+ °° 1 + N(uj 1 ) - n(w 2 ) 


/ i-oo 

doj\doj 2 - 

-OO 


UJ — — UJ2 


-^Im (-^Im <(/ 0c r|/oo->>^ 2 ) , 

where the following notation were used: 

S 0 + = /oVoL ^ = /oVot- 


(46) 


(47) 
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It is possible now to rewrite (1461) in a more convenient way 

/ L)J ^ 1 

dJ (cot -^ + tan - c*/)so<t(w')) • (48) 

The equations (1501) and (1081) constitute a self-consistent system of equations for the single-particle 
Green function of SIAM. Note that spin-up and spin-down electrons are correlated when they 
occupy the impurity level. So, this really improves the standard mean-field theory in which just 
these correlations were missed. The role of electron-electron correlation becomes much more 
crucial for the case of strong correlation. 


5 SIAM. Strong Correlation 


The simplest relevant algebra of the operators used for the description of the strong correlation 
has a similar form as for that of the Hubbard model. 28 32 79,81 Let us represent the matrix Green 
function (127j) in the following form 


G((jj) = V ( (( C ka\dlp a )) \ 

a/3 V^° Q<J l C lcr)) ((^0oo-|do/3o-))/ 

Here the operators do ar7 and are 

dicta — IT'i—cr^iai — i), 77^ — Tljcr, 77^ — (1 IT'ia'): 

^ ^ n ia ~~ 1) n ia n ia ^a/3 n iai ^ ] dj aa Oj a . 


(49) 


(50) 


The new operators di acT and dj,^ have complicated commutation rules, namely, 

[dictai ^1^0-]+ = fiijfia/3 n i—a- 


(51) 


Then we proceed by analogy with the calculations for the Hubbard model. The equation of motion 
for the auxiliary matrix Green function 


F a (u>) = 


( ((Cfcalir)) (( Cka\dl +(T )) ((cfc ff |4-<r)) \ 

{{do+a\cl a )) «d 0 +a|4+a)) «do-Hrl4-<r» 

{{do-a\cl a )) {{d 0 -a\dl +a )) ({d 0 - a \dl_ a )} 


V 




is of the following form 


EF a (uj) — I = D, 

where the following matrix notation were used 


V - £fe) -Vk 

E = I 0 (w - E 0a - U + ) 
0 0 





0 

0 ) 


II 

0 

n 0-a 

0 

II 


Vo 

0 

n 0 -J 




(52) 

(53) 

(54) 
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Here D is a higher-order Green function, with the following structure 7 11 


0 0 0 

D(uj) = ( D 21 D 22 D 23 

,D 3 i D 32 D 33 , 


( 55 ) 


In accordance with the general method of irreducible Green functions, we define the matrix irre¬ 
ducible Green function: 


= o-E(a-)) ( g “ + ° r) ' 


Here the notation were used: 


A = 


._)_ + _ ((fo-o-Cp-cr + Cp_ fr /o- C r)(nO(7 -a)) 

~ : 

((/d_ CT c p—CT + Cp_ CT /o—cr) (1 + ^0(T — n 0—a)) 

(1 - n 0 _ CT ) 

A~ + = A ++ , A + ~ = -A— 


The generalized mean-field Green function is defined by 

EF°(u)-I = 0; = 

a(3 

From the last definition we find that 

pt \\0 __ {np-a) _ ^ , Y.p V p A 


((foalfDZ = 


+ 


u-E 0 a -U + - Z P V P A ++^ + w - E 0a - C/_ 

1 — {^0—a) 


) 


)• 


w - e 0(T -U--J 2 P v p a v w - f; 0o - - u + 

{{c ka \cl)Z = {u;-e k ~\V k \ 2 F a \u))-\ 


where 


F at = 


1 - (n 0 _ ff ) 

CU F/Q(t L+ u: Eptj 17— 


(no-c 


+ 


(56) 


(57) 

(58) 

(59) 

(60) 


(61) 

(62) 


(63) 


For V p = 0, we obtain, from solution (1611) . the atomic solution F at . The conduction electron Gr< 


een 


function (J62I) also gives a correct expression for V k = 0. 


6 IGF Method and Interpolation Solution of SIAM 

To show explicitly the flexibility of the IGF method, we consider a more extended new algebra 
of operators from which the relevant matrix Green function should be constructed to make the 
connection with the interpolation solution of the Anderson model! 11 Our approach was stimulated 
by the works by J. Hubbard. 29 30 31 

Let us consider the following equation of motion in the matrix form 


k )Gpa{co) = / + Xd VpD p {u), 


(64) 


11 












where G is the initial 4x4 matrix Green function and D is the higher-order Green function: 


(G n 

G\2 

G 13 

G u \ 

G21 

G22 

G23 

G 2 4 

G31 

G32 

G33 

G34 

VC 41 

G42 

G43 

G44 ) 


Here the following notation were used 

Gu = ((Cka |cL»; G 12 = (Mil))-, 

G \3 ( (Cfccr |/oa-^O—a - ))) G\ 4 . cr)), 

G21 = «/oa|cL»; G 22 = ((Ml))-, 

G23 = ((/ocrl/oo-^O-0-)); G 2 4 = ((foa\cl a n 0 -a))', 
G31 = ((foanO- a \c\ a ))\ G32 = ((/oCT^O-crl/ocr)); 
G 33 ((/o<7^0—cr| foa^O—cr)) i G34 ((/oct^- 0 —tr))i 

^41 ((c/jo-TlO—(j|c^, cr )), Cr4 2 ((Cfc CT no-(r|/o CT ))) 

f ?43 — ((Qccr^O—o'l/ocr^'O—ex)): G *44 — ( (Cfccx^O —a | ^fco-^O—o')) • 


( 65 ) 


( 66 ) 


We avoid to write down explicitly the relevant 16 Green functions, of which the matrix Green 
function D consist, for the brevity. For our aims, it is enough to proceed forth in the following 
way. 

The equation ()64|i results from the first-time differentiation of the Green function G and is a 
starting point for the IGF approach. Let us introduce the irreducible part for the higher-order 
Green function D in the following way 

bf ] = Dp («» P) = (!. 2 > 3 ’ 4 ) ( 67 ) 

a. 

and define the generalized mean field Green function according to 


Y J F(p,k)Gpf(u) = i- ( 68 ) 

p 

Then, we are able to write down explicitly the Dyson equation (1231) and the exact expression for 
the self-energy M in the matrix form: 


M ka (u) = r 1 J2 v P v q 


P,Q 


/0 

0 

0 

Vo 


0 

0 

M 33 

M43 


0 \ 

0 

Af 34 
M44 J 


(69) 


Here the matrix I is given by 


/ 


I = 


1 
0 
0 

\(no-c 

and the matrix elements of M are of the form 


0 

1 

(no-c 

0 


0 

(no-o 
(n 0 -0 
0 


(n 0 -a)\ 
0 
0 

(no—cr)) 


(70) 


M 33 = (((4 r \p)\B[ ir \q))))W, M 34 = (((4 r \p)\Bi ir \k,q))))M 

M 43 = (((4 r \k;,p)\B[ ir \q))))W, M 44 = (((4 <r) (fc,p)|flf : 'fag))))®. 


(71) 
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Here 


A\{p) — { c p—crfoafo—cr Cp— ct/o— c/ocr) > 

2^2 (k,p) {c-ka fQ—fjCp—cr ^kaCp_ a fo—a)^ 

Bl(p) = (flcl-Jo-, ~ /ir/o-trCp-ff); ( 72 ) 

B 2 (k,p) = (4^4-cr/o-o- - C lo-/o-o- C P-o-)' 


Since the self-energy M describes the processes of inelastic scattering of electrons (c — c , / — /, 
and c — / types), its approximate representation would be defined by the nature of physical 
assumptions about this scattering. 

To get an idea about the functional structure of our generalized mean field solution (f 68 l) . let us 
write down the matrix element G^ F : 


G 33 F = ((foan 0 -a\foa n O-a)) = 
(no -a) 


w - ef F - U - S mf (uj) - Y(u) 
(n 0 -a)Z(uj) 


+ 


Z(u) = S(u )^2 


(w - ef F - U - S MF (u) - Y(u))(u - E 0a - 5(w)) ’ 

UZ{u) 

LO - E 0a ~ S(U))' 

32 

, V P^ 

P 


Y(u) = 


Vpha +E^^ + s(^"+E^ 


OJ - 6 
p P 


MF 


co - e 1 ' 
p p 


(73) 

(74) 

(75) 


Here the coefficients L 41 , L 42 , L 31 , and L 32 are certain complicated averages (see definition (|67[) ) 
from which the functional of the generalized mean field is build. To clarify the functional structure 
of the obtained solution, let us consider our first equation of motion (16411 . before introducing the 
irreducible Green functions (1671) . Let us put in this equation the higher-order Green function 
D = 0. To distinguish this simplest equation from the generalized mean field one (1681) . we write 
it in the following form 

^F(P,k)G°(p,u J ) = i. (76) 

p 

The corresponding matrix elements which we are interested in here read 


1 - (riQ-cr) 


+ 


G° 22 = ((foolfl)) = 

(no-a) 


0 i-Eoa-S{uj) u-Eov- S(u)-U' 


G 33 = ((foan 0 -a\fo a n O-a-)) = 


(n 0 -o 


u-E 0(T - S(u)-U’ 

<4 = ((fo.n 0 -Ml)) = G§3- 


(77) 


(78) 

(79) 


The conclusion is rather evident. The simplest interpolation solution follows from our matrix 
Green function (|65l) in the lowest order in V, even before introduction of generalized mean field 
corrections, not speaking about the self-energy corrections. The two Green functions G® 2 and G 33 
are equal only in the lowest order in V. ft is quite clear that our full solution (|241) that includes 
the self-energy corrections is much more richer. 

ft is worthwhile to stress that our 4x4 matrix generalized mean field Green function (|65l) gives 
only approximate description of suitable mean fields. If we consider more extended algebra of 
relevant operators, we get the more correct structure of the relevant generalized mean field. 
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7 Quasiparticle Dynamics of SIAM 


To demonstrate more clearly the advantages of the irreducible Green functions method for SIAM, 
it is worthwhile to emphasize a few important points about the approach based on the equations- 
of-nrotion for the Green functions. To give a more instructive discussion, let us consider the 
single-particle Green function of localized electrons G a = ((/oer|/(L-))- The simplest approximate 
’’interpolation” solution of SIAM is of the form: 

1 U(n 0 -a) 


G a (oj) = 


u-E 0a - S(w) + (w - E 0a - S(w) - U)(u - E 0rT - 5(w)) 

1 {IT'O—a ) (^ 0 —it) 


+ 


u-E 0a - S(w) u-E 0a - S(lu)-U' 

The values of n a are determined through the self-consistency equation 


n a = (n 0r j) = -~ dE n(E)lm G C7 (E, n a ). 

7 T J 


The atomic-like interpolation solution 

g ^ ) = EE^L + 

W ~ E Oo- 

GJu>) = 


reproduces correctly the two limits: 
(no-a) 


u-E 0rT -U' 

1 


where 


w - E 0c7 - S(u) ’ 

\V k \ 2 


for V = 0, 
for U = 0, 


s ( w) = 


to -e k 


(80) 

(81) 

(82) 

(83) 

(84) 


The important point about equations (1821) . (1831) is that any approximate solution of SIAM should 
be consistent with it. Let us remind how to get solution (182p . It follows from the system of 
equations for small-U limit: 

(a; - E 0a - 5'(cj))((/ 0CT |/ 0 t (T ))a; = 1 + U ((/ OCT n 0 - CT |/o CT ))a;, 

(u - E 0a - [/)((/o <T no- ff |/ 0 t (7 ))c ~ 

(^0—cr) T Lfc((cfco-^'0a-|/d a -))a;; 


k 

{u - e fc )((c fc(T no- (J |/ 0 t (T ))a; = ^((/oo-no-o-l/do-))^ 

Note that the equations ([551) and (1551) are approximate; they include two more terms. 
We now proceed further. The starting point is the system of equations: 

(u - E 0a - S{uj)){(f 0rT \foa)) = 1 + U{(f 0a n 0 -a\foa))’ 

(u - E 0(7 - C/)((/ 0o -n 0 - t r|/d (7 )) = (n 0 -a) + ^2 V k (((Cfc CT ™ 0 - CT |/(L)) ~ 


{{ck-afl-akAfla)) T ((Cfc-o-Zocr/o-trl/dtr))^ • 


(85) 

( 86 ) 

(87) 

( 88 ) 


Using a relatively simple decoupling procedure for a higher-order equation of motion, a qualita¬ 
tively correct low-temperature spectral intensity can be calculated. The final expression for G for 
finite U is of the form 

1 

I Oo-/' 


((ML)) = 


+ 


uj-E 0(T -S{u;) + US 1 {uj) 
U(no-a) + UFi(uj) 

K{uj) (uj - E 0(7 - S(u) + C75i(w)) : 


(89) 
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where F\, Si, and K are certain complicated expressions. We write down explicitly the infinite U 
approximate Green function: 


((foMl)) = 


1 (^0—<x) Frr (^) 

u; - E 0a - S(u) - Zl(u)' 


The following notation were used: 


v 2 T, 


q,k 


( C l-qCk-a) 
u - e k 


F a 




k 


(fo~a C k-o-) 
u - e k 


k 


(fo-cr^k-a ) 

uj- e k 


(90) 


(91) 

(92) 


We putted here V k — V for brevity. The functional structure of the single-particle Green function 
(1891) is quite transparent. The expression in the numerator of (1891) plays the role of an effective 
dynamical mean field, proportional to (f})_ a Ck-a) ■ In the denominator, instead of bare shift S{u) 
AMD we have an effective shift S 1 = S(lj) + Z\{u). The choice of the specific procedure of 
decoupling for the higher-order equation of motion specifies the selected generalized mean fields 
and effective shifts. 


8 Complex Expansion for a Propagator 

We now proceed with analytic many-body consideration. One can attempt to consider a suitable 
solution for the SIAM starting from the following exact relation derived in paper : 11 

((fo*\fl)) = g 0 + 9°Pg 0 , (93) 

g° = (u J -E 0 a-S(u))- 1 , (94) 

P = U (n 0 -a) + U 2 ((f 0a n 0 -a\fln 0 -a)). (95) 

The advantage of the equation (fMl) is that it is a pure identity and does not includes any 
approximation. If we insert our generalized mean field solution (1901) into (|93D . we get an essentially 
new dynamic solution of SIAM constructed on the basis of the complex (combined) expansion of 
the propagator in both U and V parameters and reproducing exact solutions of SIAM for V = 0 
and U = 0. It generalizes (even on the mean-field level) the known approximate solutions of the 
Anderson model. 

Having emphasized the importance of the role of equation (1931) , let us see now what is the best 
possible fit for higher-order Green function in (1951) . We consider the equation of motion for it: 

(U - E 0 a ~ U)({f 0 an 0 -a\foa n O-<r)) = (™ 0 -<r) + (96) 

^ ^ kfc(((Cfco-^0 —o'I/qo-^'O— g)) T 
k 

((4-<r/o«T/o-ff|/o <T nO- ( r)) - ((Ck-afo^afoa\foa n O-a)))- 

We may think of it as defining new kinds of elastic and inelastic scattering processes that con¬ 
tribute to the formation of generalized mean fields and self-energy (damping) corrections. The 
construction of suitable mean fields can be quite nontrivial, and to describe these contributions 
self-consistently, let us consider the equations of motion for higher-order Green functions in the 
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r.h.s. of ([MD 


(a; - efc)((cfc 0 -n 0 _ (T |/Q (T no_ 0 -)) = 

v {{foano-alfl^riQ-o)) + 

^ ] V((( c kcrfo—a c p—Q-\foa n O—a)} ~~ ( ( c ka c p— cr/o —a \ /cL n 0— a)) )j 

V 

^k ^Ocr “ 1 “ ^ 0 —cr) ((Cfc—o'/ , Q_ (T /ocr|/oa-^O—( t)) 

{fo—a^k—a^Oa) 

v {{foano-alfl^no-o)) + 

^VdiCk-rf^Cp^f^no-a)) ~ ((Ck^cl^foalfoa^-a))), 

V 


(97) 


(98) 


(w + e k - E 0rT - Eo-v - f7)((4_ CT / OCT /o—al/o^no— ff }) (99) 

= -( cl_ a foafLfo-a) + 

V{{foano-a\foa n O-cr)) + 

^2V(((cl_ a C pa fo- a \f^lO-a)) + ( Ocr Cp—cr | fQu^O-cr )) ) ■ 

P 


Now let us see how to proceed further to get a suitable functional structure of the relevant solution. 
The intrinsic nature of the system of the equations of motion (f97l) - (f99l) suggests to consider the 
following approximations: 

(w - e fc )((c fco .no_ CT |4 0 .no_ 0 -)) « V((f 0a n 0 - a \f^ a n 0 - a }), (100) 

(^ ffc Eo a + Eo- a )((ck-afo- a focr\foa n O-'?)} ~ (fo-cr^k—a^Oa) 

-Vdifoano-alf^no-a)) - ((c fc _ ff 4- CT /oa|/(L™o-a)}), (101) 

(w + e k - E 0a - Eq— fy - U)((cl_' T fo a fo-< T \fL n o-'T)) ~ -{c\_ a foafl a k-o) + 

V (((/oct^O-ctI/oct^O-ct)) + {{cl- a foaC k -cr\fL n O-<7)))- (102) 


It is transparent that the construction of approximations (HOOD - (11021) is related with the small -V 
expansion and is not unique, but very natural. As a result, we find the explicit expression for 
Green function in (1951) 


((/oct^O— o -1 /oct^O— a)) 


Here the following notation were used: 


{no- a ) - Fl(u) 
lj - E 0a - U - Siiu)' 


(103) 


Si(w) = S(w) 

+ VlH| 2 (_-_|_ - _) 

■“ w — f-k ~ -E'oo- + T^o-o- w + efc — Eo a — Eo_ a — U 

El = ]T(HF 2 + V 2 F 3 ), 

k 

p _ i^k—crfocr f 0(j/o —<t ) + (fo— a Ck—cr n Oa) 

(d Eqq- Eq—q- U UJ Eqq- Eq—q- 

((c fe _ CT c{._ CT /o CT |/ 0 t fT no— a)) {(c\_ a fo (7 C k -a\fo CT n O-a)) 

r 3 = --- 1 ----- . 

oj a k Eo a + E§— a uj + € k Eo a -Eq—< j U 


(104) 


(105) 

(106) 
(107) 
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Now one can substitute the Green function in (1951) by the expression (11031) . This gives a new 
approximate dynamic solution of SIAM where the complex expansion both in U and V was 
incorporated. The important observation is that this new solution satisfies both the limits (I82|) . 
For example, if we wish to get a lowest order approximation up to U 2 and V 2 , it is very easy to 
notice that for V = 0: 


((fo<j c k—a C k—a l/ocr^O— a)) 

{{Ck-oc\_ a ha\fl a n 0 -a)) 


( C l_c r c k—cr)( n 0—cr) 

W ~ Eoa ~ U 
{c k - a c\_ a )(n 0 - a ) 
W ~~ Eq<t ~ U 


(108) 

(109) 


This results in the possibility to find explicitly all necessary quantities and, thus, to solve the 
problem in a self-consistent way. The main results of our IGF study is the exact Dyson equation for 
the full matrix Green function and a new derivation of the generalized mean field Green functions. 
The approximate explicit calculations of inelastic self-energy corrections are quite straightforward 
but tedious and too extended for their description. Here we want to emphasize an essentially new 
point of view on the derivation of the generalized mean fields for SIAM when we are interested in 
the interpolation finite temperature solution for the single-particle propagator. Our final solutions 
have the correct functional structure and differ essentially from previous solutions. 

In summary, we presented here a consistent many-body approach to analytic dynamic solution 
of SIAM at finite temperatures and for a broad interval of the values of the model parameters. 
We used the exact result (1931) to connect the single-particle Green function with higher-order 
Green function to obtain a complex combined expansion in terms of U and V for the propagator. 
We reformulated also the problem of searches for an appropriate many-body dynamic solution 
for SIAM in a way that provides us with an effective and workable scheme for constructing of 
advanced analytic approximate solutions for the single-particle Green functions on the level of the 
higher-order Green functions in a rather systematic self-consistent way. This procedure has the 
advantage that it systematically uses the principle of interpolation solution within the equation- 
of-motion approach for Green functions. The leading principle, which we used here was to look 
more carefully for the intrinsic functional structure of the required relevant solution and then to 
formulate approximations for the higher-order Green functions in accordance with this structure. 
Of course, there are important criteria to be met (mainly numerically) , such as the question 
left open, whether the present approximation satisfies the Friedel sum rule (this question was 
left open in many other approximate solutions). A quantitative numerical comparison of self- 
consistent results e.g. the width and shape of the Kondo resonance in the near-integer regime of 
the SIAM would be crucial too. In the present consideration, we concentrated on the problem of 
correct functional structure of the single-particle Green function itself. 


9 The Improved Interpolative Treatment of SIAM 

For better understanding of the correct functional structure of the single-particle Green function 
the development of improved and reliable approximation schemes is still justified and necessary, 
and an effective interpolating approximations are desirable. The present section is devoted to the 
development of an improved interpolating approximation 10 11 for the dynamical properties of the 
SIAM. We will show that a self-consistent approximation can be formulated which reproduces all 
relevant exactly solvable limits of the model and interpolates between the strong- and the weak- 
coupling limit. This approach is complementary to the one described above. 

We start by considering the equations of motion for the Fourier transformed Green function, 

POO 

G a {u) = {{fo\ft))u> = -i dt ex.p(iujt)([foa{t), /ocr]+) : (HO) 

J 0 
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(u-E a - 5(w))«/ ff |/t)) u = 1 + U((f a n- a \ft)) u = 1 + X<,{u,)((U\fi))u 


( 111 ) 


Here the the quantity E a (uj) may be conditionally interpreted as the one-particle self-energy and 

IV? 


S( w) = X] 


W - e k 


( 112 ) 


We want to develop an interpolating solution for the SIAM, i.e. a solution which is applicable in 
both, the weak-coupling limit (and thus the exactly solvable band limit) and the strong-coupling 
limit (and thus the atomic limit). As it was shown earlier, the simplest approximative interpolating 
solution has the form: 


, + 


(n-a) 


uj — E a — S (oj) uj — E a — S(uo) — U 


(113) 


Here {n- a ) denotes the occupation number of /—electrons with spin a. This is just the analogue 
of the Hubbard III approximation 28 for the SIAM. As for the Hubbard model, however, Fermi 
liquid properties and the Friedel sum rule, which hold for the SIAM at least order by order within 
the [/-perturbation theory, are violated within this simple approximation. 

An approximation, which automatically fulfills Fermi liquid properties and sum rules, is provided 
by the self-consistent second order [/-perturbation treatment (SOPT) and is given by: 

E a (iujn) = U(n- a ) - f^r j G a (iu n + iv)G- a {iui - iv)G _ ff (wi). (114) 

Here wi(z/) denote odd (even) Matsubara frequencies and (3 = 1 /ksT. One of our goals is to find 
some way to incorporate this SOPT into an interpolating dynamical solution of the SIAM. This 
means that the approximation for the self-energy shall be correct up to order U 2 perturbationally 
around the band limit U = 0 and also the atomic limit V = 0 shall be fulfilled. This is the 
case for the SOPT around the Hartree-Fock solution, but only for the symmetric SIAM. For the 
general situation (position of the Fermi level relative to E a and E a + U) a heuristic semi-empirical 
approach only for constructing such an approximation has been discussed in literature Here our 
intention is to take into account the self-consistent-SOPT. Furthermore, the approximation shall 
not only fulfill the atomic limit V = 0, but it shall be correct up to order V 2 in a strong-coupling 
expansion around the atomic limit. 

The self-consistent inclusion of contributions in second (and fourth) order perturbation theory 
around the atomic limit is, in particular, important to properly account for the Hondo effect 
within the SIAM (Hondo temperature scale) and to reproduce the correct antiferromagnetic be¬ 
havior in the strong-coupling limit of the Hubbard model. Especially the calculation of some 
magnetic properties for the Hubbard model and the well known Hondo effect for the SIAM shows 
the importance of second (and fourth) order perturbation theory around the atomic limit. 

It was already mentioned that during the last decades several different refined many-body tech¬ 
niques have been applied to the SIAM, and many of these approaches are strong-coupling treat¬ 
ments around the atomic limit and can be classified as being correct up to a certain power in the 
hybridization V. When applied to the calculation of static properties many of these treatments, 
give reasonable results. But for the many-body dynamics the results of most of these approxima¬ 
tions are not fully satisfactory, in particular as Fermi liquid properties and sum rules are violated. 
Furthermore, when applied to the finite-// SIAM none of these approximation schemes reproduce 
the SOPT, i.e. these approaches are not correct in the weak-coupling limit up to order U 2 . 

To construct the interpolating approximation 7 10 11 8 for the SIAM fulfilling all desired proper¬ 
ties mentioned above we start from the equation of motion for the higher order Green function 
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( 115 ) 


{{fan- a \fl))u- 

(u-E a - S(u) - U)({fan-a\fl)) u = (n~a) - 


With 

[G^w )] -1 = u-E <r - S(u) 

and the self-consistent summation 


[G^°^(a;)] 1 G a ( uj) — 1 + E a (uj)G a (uj), 


we derive from this equations of motion the following exact relation 

’ [U 1-(U -Y,„(u))G,(li) 


Here the definition 


{{faf-<Tf-a\fin- a ))uj 


-Z{u) 


G a (u>) 

1 + Tt a (ui)G a -(ui) 


was introduced. 

Applying the equations of motion to the higher-order Green function 


({faf-crf-a\fa n -a})u 


(116) 


(117) 


(118) 


(119) 


one obtains for the function Z(z) the exact equation 


Z(u) = V V] { Gi a (k) - G 2 CT (k) + —— [G 3 CT (k) - G 4 rJ ( k)]l, (120) 

^ l u ~ e k J 

with k = (k, uj) and 

Glaik) = (121) 

G2*( k) = (ifA-J-* \ftn-°)U (122) 

G 3 c r(k) = ^((c k0 -/- 0 -4_ CT |/ln_ (T )) w , (123) 

q 

GAa (k) = (124) 

q 

Self-consistency in the perturbation theory defines the Green function: 

[G^lcu)] 1 G a {uj) = 1 + S 0 .(cu)G fr (o;), (125) 


and leads to an infinite order resummation resulting in a self-consistent approximation. 

In general, there are several possibilities to incorporate self-consistency, but most of these possi¬ 
bilities lead once more to an approximation being exact up to order V 2 but not reproducing the 
weak-coupling limit. To be exact up to order V 2 it is justified to replace the higher order Green 
functions on the right hand side of Eq. (112011 by their lowest order contributions, which are given 
by 


Gla (k) 

G 2 a( k) 
GW(k) 


V 


e k - E_ a - U 


{n-a ) L/k - f(E-a + u)] + (n-a) [1 - / k ] (n-o) [1 - / k ] 


CU 6 k Efy + E—Q 


V 


E _ 0 


(1 - (n-o))[f k - f (E^ a )\ + [1 - /k](n_ CT ) 
w + e k — E cT — E_ a — U 


uj — E a — U 
(n-o)[ 1 - /k] 


+ 0(H 3 ), 


uj — E a — U 


+ 0{V% 


OiV 2 ), G 4 a (k) = 0(V 2 ), 


(126) 
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leading to a finite order V 2 perturbation expansion of the self-energy (11171) . Here f{E) = {exp[(E— 
n)/ksT] + 1 } _1 is the Fermi function, // the chemical potential and /k = /(e k). 

For the higher order Green functions Gi a { k) (i = 1,... ,4) one can find an approximation which 
reproduces the exact relations (11261) in lowest order in V and is simultaneously exact in lowest 
order in U (when Wick’s theorem is applicable). One possibility for such an approximation is 
given by: 


G i ff (k) = 


x 

G 2ct (k) = 


x 


X 


7 -rr-r f<y\ n -°fl))iun+iv(( c k-<y\naf-a)) 

(n- a )(n a n- a ) 


— CT//IUJ1—IV 






~P 


-2 


(1 


(U— a ) 

EG—rr— f-°S-Jl)) 




{{f-c\n<rf'L a ))iu 1 -i v - «/ ff |n-jJ))i Wn+l G(/-G/jJ/I (J ))»u, 1 -u 




(/o-/o-cj c _ 0 ./_< 7 ) u £ ___ | j.f.. 


1 - {n_ a ) 


■(( Un~ a \fl))i 


(127) 


(128) 


and the Green function are decoupled according to the theorem of Wick. Since the 

approximation does not violate the theorem of Wick for small U, it automatically satisfies the 
SOPT, i.e., expanding Eq. (I117D for small U up to second order in U leads to the SOPT for 
the self-energy. Also the V 2 - limit is not violated since the Green function G^a^G^ are itselves 
proportional to V 2 , leading in Eq. (11201) to V 4 terms. Therefore, our approximation leads to an 
expression for the self-energy of the SIAM, which is exact at least up to order U 2 in a weak cou¬ 
pling expansion and up to order V 2 in a strong coupling expansion. The structure of the chosen 
approximation (I127D and (J128D and of the decoupling for the Green function G ^ a , G 4 a according 
to the theorem of Wick has a similar analytical structure as the SOPT, (which can be calculated 
numerically very fast and accurate). Hence the explicit numerical calculations within this treat¬ 
ment are of the same order of complexity as those of the self-consistent-SOPT calculations. 
Notice that in principle it is possible to systematically improve the above approximation. Since 
the self-consistent summation (11171) . (11201) is formally exact, the next step would be the similar 
construction of an approximation for the Green functions G^a^G^ (and for Green functions of a 
similar structure occurring in a further application of the equations of motion to the Green func¬ 
tions G\ a ,G 2 a) being exact in order V 2 and simultaneously satisfying the theorem of Wick; as 
the Green functions G^a etc. have already a prefactor V 2 in (I120D this leads to an approximation 
for S and thus the self-energy E CT (u;) being exact up to order V 4 in the strong-coupling limit and 
simultaneously in order U 2 in the weak-coupling limit. Furthermore, already from the structure of 
the exact equation (11171) it is clear that our new approximation can be considered as a systematic 
improvement of the Hubbard-III approximation (I113D . which is known to be reasonable concerning 
the high-frequency behavior of the dynamical quantities and concerning the reproduction of the 
metal-insulator transition in the Hubbard model. 

The improved approach goes beyond the Hubbard-III approximation 28 including all self-energy 
contributions in order U 2 and thus reproducing the SOPT. This is important to fulfill the Fermi 
liquid properties at least for small 17, and in this respect the approach should be as good as the 
related attempts. 

Ou the other hand, the new approach is also exact up to order V 2 and is, therefore, as good as 
standard equations of motion decouling procedures are, which qualitatively describe important 
items like Kondo peak, Kondo temperature scale, etc. 
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When interpreting these standard equations of motion decouplings as generalized mean-field treat¬ 
ments, because the decoupling consists in a replacement of a higher order Green function by a 
product of an expectation value with a lower order Green function, our new approximation can be 
considered to be a kind of dynamical mean-field approximation, because the approximation (1127(1 . 
(I128D consists in the replacement of a higher order Green function by combinations of products 
of (time-dependent) lower order Green functions. 

Finally the approach is not a completely uncontrolled approximation, as it is exact up to certain 
orders (V 2 , U 2 ) of systematic perturbation theory. It is, however, as any self-consistent approx¬ 
imate treatment is, uncontrolled in the way it takes into account infinite order resunmiations of 
arbitrary order in U and V by the self-consistent requirement, which is unavoidable to reproduce 
both limits. 

In summary, an improved interpolating approximation for the SIAM has been developed, which 
recovers the exactly solvable limits V = 0 and U = 0 and which is even more at least correct up 
to order V 2 in a strong-coupling expansion and simultaneously up to order U 2 in a weak-coupling 
expansion. 


10 Quasiparticle Many-Body Dynamics of PAM 

The main drawback of the Hartree-Fock type solution of PAM ([6|) is that it ignores the correlations 
of the ””up” and ’’down” electrons. In this section, we will take into account the latter correlations 
in a self-consistent way using the IGF method. We consider the relevant matrix Green function 
of the form (cf. (f27l) ) 


Q( 0 j\ — ( (( c fco’l c feo-)) ({ c ka\ft a )) 

V«/*dcL» «iWI/b>, 

The equation of motion for Green function (11291) reads 


(129) 


(W - € k ) —14 

— Vfc (U! — E k ) 


(( Cka\c\ a )) (( C k(r \fL » 
{{fka\c\. a }} {{fka\flo)) . 



0 0 
J w ///iid 


«AcL» (WU) 


(130) 


where A = f k+pa f j )+q - a fq-a- According to IGF method the definition of the irreducible parts in 
the equation of motion (11301) are given by 

{lr) {{fk+ P afl +q -afq-'j\ c \o)) = ((fk+pa fl+q-a /<?-d c L» “ 5 pfi{ n q-<T){{fka\c\ a )), (131) 

M {{fk+pafl +q -afq-a\fL)) = {{fk+pafl +q - a fq-a\fl a )) ~ Sp, 0 {n q - a ) ((f ka |/^)). (132) 

After substituting these definitions into equation (11301) . we obtain 

(“ - a) -Vk \ (((Ckv l4„>> <<Cto_ 

-u (c-E„w)H((/ to |cL)> «/toi/L))r 


1 0 \ TTAT~^-\~^ ( 0 0 


0 1 


(133) 


In the the following the notation will be used for brevity 


E a (k) = E k - Un - a ; n- a = {fl.fk-a)- 


(134) 
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The definition of the generalized mean field Green function (which, for the weak Coulomb cor¬ 
relation U, coincides with the Hartree-Fock mean field ) is evident. All inelastic renormalization 
terms are now related to the last term in the equation of motion (|133f) . All elastic scattering ( or 
mean field) renormalization terms are included into the following mean-field Green function 


/(w-e t ) -V t \{((ctc l<b»° (tetol/L»°) _ A °V 

V -U (* - B.(k))) («/ to i4„»° «/ to i/L)>V v» V' 


It is easy to find that (c/. (1321) and (|33j)l 



<</^i/L>>° = ( 


(135) 

V u) - e k ) 


<fe.ii»° = 

( . . I^l 2 \- X 

V 6k uj — E a (k)) ' 

(136) 


At this point, it is worthwhile to emphasize a significant difference between both the models, 
PAM and SIAM. The corresponding SIAM equation for generalized mean field Green function 
(13TD reads 



-V p 5 pk \ ({{cka\c ] ko ))° 

— Efta — Un-a)J ^((/oo-lcj^)) 0 


((f0a\fl))°J 



(137) 


This matrix notation for SIAM shows a fundamental distinction between SIAM and PAM. For 
SIAM, we have a different number of states for a strongly localized level and the conduction 
electron subsystem: the conduction band contains 2 N states, whereas the localized (s-type) level 
contains only two. The comparison of (11371) and (11351) shows clearly that this difficulty does not 
exist for PAM : the number of states both in the localized and itinerant subsystems are the same, 
i.e. 2N. 


This important difference between SIAM and PAM appears also when we calculate inelastic scat¬ 
tering or self-energy corrections. By analogy with the Hubbard modelj 32l ^#Sll the equation of 


motion (I133D for PAM can be transformed exactly to the scattering equation of the form 
Then, we are able to write down explicitly the Dyson equation (1231) and the exact expression for 
the self-energy M in the matrix form: 




0 0 \ 
0 M 22 ) 


Here the matrix element M 22 is of the form 


M22 — M ka (u) — 


El 

N 2 


^( { " ) ((/ fc+pCT / P Va/' ? --l/r- CT /r+- CT / fe t +SCT )) (ir) ) (P) - 

pqrs 


(138) 


(139) 


To calculate the self-energy operator (I139D in a self-consistent way, we proceed by analogy with 
the Hubbard model. Then we find both the expressions for the self-energy operator 32 79 81 by 
iteration procedure. 
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11 Quasiparticle Many-Body Dynamics of TIAM 


Let us see now how to apply the results of the preceding Sections for the case of TIAM Hamiltonian 
(|7|) . The initial intention of Alexander and Anderson 3 was to extend the theory of localized 
magnetic states of solute atoms in metals to the case of a pair of neighboring magnetic atoms . 8 40 
It was found that the simplified model based on the idea that the important interaction is the 
diagonal exchange integral in the localized state, which is exactly soluble in Hartree-Fock theory 
for isolated ions, is still soluble, and the solutions show both ferromagnetic and antiferromagnetic 
exchange mechanisms. 

Contrary to that, our approach go beyond the Hartree-Fock approximation and permits one to 
describe the quasiparticle many-body dynamics of TIAM in a self-consistent way. 

We again consider the relevant matrix Green function of the form (cf. (|27|) ) 


f G\\ G 12 G 13 > 

G(lo) = | G 21 G 22 G 23 | = 

^31 G32 G33, 


«cUi» (Mfl)) «cU/L»\ 

,t u // r. 1 rt \\ // 1, 1 f t 

k 

,t 


((/idcL)) ((hMD) ((fi*\f&) 


K ((f2a\4J) ((f2a\fl)) ((f2a\fl))J 

The equation of motion for Green function (11401) reads 

E 

p \ -V2 P 

0 0 N 

0 1 Ol+Cp^lcU) {(^i|/L)> ((AilfD) 

0 o V((a 2 |cL» «^2i/L» ((^2i/k», 


/ (w €p ) dp^ 

1 Ipfipk 

Vlpfipk \ 

( Gu 

G12 

G i 3 \ 

-Vip 

7 v( W B/Qu) 

-V 12 

G21 

G22 

g 23 

\ -V2p 

-V21 

Jj(i 0 — Eoa)) 

\G 3 1 

G32 

G33J 


The notation are as follows 


M — flo-fl-vfl-a, M — f 2 af\^ a f 2 -a- 


In a compact notation, the equation (I141D has the form 

J2F(p,k)G pk (u;) = i + UD p (u). 


(140) 


(141) 


(142) 


(143) 


We thus have the equation of motion (I143D which is a complete analogue of the corresponding 
equations for the SIAM and PAM. After introducing the irreducible parts by analogy with the 
equation 


{ir) ((hafl-Jl- f r\B)) u = ((flofl_Jl-a\B)) u - (n 1 _ (r )((/ 1 <T | J B)) w , 
{ir) ((f 2 afLah-a\B)) w = ((f 2 afLj 2 -a\B))^ ~ (n 2 . a ) <(/ 2(T \B)) U , 


(144) 


and performing the second-time differentiation of the higher-order Green function, and introducing 
the relevant irreducible parts, the equation of motion (I143D is rewritten in the form of Dyson 
equation (1231) . The definition of the generalized mean field Green function is as follows 
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The matrix Green function (11451) describes the mean-field solution of the TIAM Hamiltonian. The 
explicit solutions for diagonal elements of G° are 


((Cka\ cL»“ = (w - €fc - ^ ^ | Un _^ ~ , 

((fiMl))°u = (<*> - (E 0 a ~ Uri-a) - S(ui)) - A 22(k,u>)y\ 

{{f 2 o\fl))l = (w - (Eo, - U Tl-a) - S(u)) - A 33 (k,u)y\ 
Here we introduced the notation 


UJ (Etycr Un — fy'j 

UJ (Tiber Un —fj) 


UJ (-^Ocr 7777 .—o')' 

_ V 21 V 12 _' _1 

uj (E 0a Un —o - )- 


A22 (k,Uj) — (A 2 l(w) + V12) (A2I (uj) + V21) w — (Eocr — Un- a ) — 


A- 33 (k,Uj) — (A 12 (uj) + V 21 )(A 12 (w) + V 12 ) UJ — ( Eq, a — Un- a ) — 


-1 


Yh P \ v 2 pY 


UJ - e r, 


EJ^lV 1 


uj — e 


A 12 — A 21 — ^ 


p 

1 lpVbp 


uj — e 


p 


The formal solution of the Dyson equation for TIAM contains the self-energy matrix 

/0 0 0 \ 

M = Jo M 22 M‘23 I , 

\0 M 32 M33 / 


where 


(146) 

(147) 

(148) 


(149) 

(150) 

(151) 

(152) 


(153) 


M 22 = U\^((f la n 1 . tT \flrii-a)) (ir) y, 

(154) 

M 32 = U\^((f 2a n 2 .Mlni-a)) {ir) ) p , 

(155) 

M 23 = U\^\(f lani .Mln2-a)) {ir) ) p , 

(156) 

M 33 = U\W((f 2(T n 2 ^\fln 2 -*)) {ir) ) p . 

(157) 


To calculate the matrix elements (11541) . the same procedure can be used as it was done previously 
for the SIAM (1441) . As a result, we find the following explicit expressions for the self-energy matrix 
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elements fcf. (l46l) 


/ +oo 

duidu}2 

-OO 


1 + N(uji) — n(u; 2 ) 

OJ — 0J\ — UJ2 


-i/m((5r|5+)) Wl )(-i/m((/ 4 |/ 1 t i )) W2 ), 


/ +oo 

dtjJiduj• 
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/ +oo 
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-OO 


1 + N(u>i) - n(u> 2 ) 

OJ — 0J\ — UJ2 
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r+oo 

Af 23 (u;) = U 2 doj\duj2 
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1 + N(oj i) — n(cj 2 ) 
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r ^7m«5+|5 1 -)) W i) (-^m((/i t |/ 2 f t )) W2 

Here the following notation were used: 

^ = 4/h; -S f!M i = 1,2. 


( 158 ) 


(159) 


(160) 


(161) 


For A /33 we obtain the same expressions as for M 22 with the substitution of index 1 by 2. For 
'M 

A /32 we must do the same. It is possible to say that the diagonal elements M 22 and A /33 describe 
single-site inelastic scattering processes; off-diagonal elements A / 2 3 and M ^2 describe intersite in¬ 
elastic scattering processes. They are responsible for the specific features of the dynamic behavior 
of TIAM ( as well as the off-diagonal matrix elements of the Green function G°) and, more gener¬ 
ally, the cluster impurity Anderson model (CIAM). The nonlocal contributions to the total spin 
susceptibility of two well formed impurity magnetic moments at a distance R can be estimated as 


Xpair ~ «Sr|S 2 + )) ~ 2 x - 127r£ F (— 

9dB 


2 cos(2 kpR) 

(.k F R ) 3 


(162) 


In the region of interplay of the RKKY and Hondo behavior, the key point is then to connect the 
partial Hondo screening effects with the low temperature behavior of the total spin susceptibility. 
As it is known, it is quite difficult to describe such a threshold behavior analytically. However, 
progress is expected due to a better understanding of the quasiparticle many-body dynamics both 
from analytical and numerical investigations. 


12 Conclusions 

In summary, we presented in this paper in terse form a general technique how a dynamical solu¬ 
tion for SIAM and TIAM at finite temperatures and for the broad interval of the values of the 
model parameters can be constructed in the spirit of irreducible Green functions approach. We 
used an exact result to connect the single-particle Green function with the higher-order Green 
function to obtain an complex expansion in terms of U and V for the propagator. This approach 
provides a plausible yet sound understanding of how structure of the relevant dynamical solution 
may be found. Hence this approach offer a both powerful and workable technique for a systematic 
construction of the approximative dynamical solutions of SIAM, PAM and other models of the 
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strongly correlated electron systems. 

In short, the theory of the many-body quasiparticle dynamics of the Anderson- and Hubbard-type 
models at finite temperatures have been reviewed. We stressed an importance of the new exact 
identity relating the one-particle and many-particle Green functions for the single-impurity An¬ 
derson model: G = go + goPgo■ 

The application of the IGF method to the investigation of nonlocal correlations and quasiparticle 
interactions in Anderson models 8 has a particular interest for studying of the inter-site correlation 
effects in the concentrated Kondo system and other problems of solid state physics . 40 65 A compar¬ 
ative study of real many-body dynamics of single-impurity, two-inrpurity, and periodic Anderson 
model, especially for strong but finite Coulomb correlation, when perturbation expansion in U 
does not work, is of importance for the characterization of the true quasiparticle excitations and 
the role of magnetic correlations. It was shown that the physics of two-impurity Anderson model 
can be understood in terms of competition between itinerant motion of carriers and magnetic 
correlations of the RKKY nature. This issue is still very controversial and the additional efforts 
must be applied in this field. 

The many-body quasiparticle dynamics of the single-impurity Anderson Model was investigated 
by means of the equations of motion for the higher-order Green functions. It was shown that 
an interpolating approximation, which simultaneously reproduces the weak-coupling limit up to 
second order in the interaction strength U and the strong coupling limit up to second order in 
the hybridization V (and thus also fulfills the atomic limit) may be formulated self-consistently. 
Hence, a new advanced many-body dynamical solution for SIAM has been developed, which re¬ 
covers the exactly solvable limits V = 0 and U = 0 and which is even more at least correct up to 
order V 2 in a strong-coupling expansion and simultaneously up to order U 2 in a weak-coupling 
expansion. 

Further applications and development of the technique of the equations of motion for the Green 
functions were described in Refs . 34 86,8 ' 88 89 90 91 HIOHIE SI These applications illustrate some of 
subtle details of this approach and exhibit the physical significance and operational ability of the 
Green function technique in a representative form. 

This line of consideration is very promising for developing the complete and self-contained theory 

main results reveal the fundamental importance of the adequate definition of generalized mean 
fields at finite temperatures, that results in a deeper insight into the nature of quasiparticle states 
of the correlated lattice fermions and spins. We believe that our approach offers a new way for 
systematic constructions of the approximate dynamic solutions of the Hubbard, SIAM, TIAM, 
PAM, spin-fermion, and other models of the strongly correlated electron systems on a lattice. 
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